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Abstract. The two independent "plus" and "cross" polarization waveforms associated 
with the gravitational waves emitted by inspiralling, non-spinning, compact binaries 
are presented, ready for use in the data analysis of signals received by future 
laser interferometer gravitational- wave detectors such as LIGO and VIRGO. The 
computation is based on a recently derived expression of the gravitational field at 
the second-post-Newtonian approximation of general relativity beyond the dominant 
(Newtonian) quadrupolar field. The use of these theoretical waveforms to make 
measurements of astrophysical parameters and to test the nature of relativistic gravity 
is discussed. 



1. Introduction 

Two large-scale laser interferometer detectors of gravitational waves are now under 
construction in the U.S. (LIGO experiment [1]) and in Europe (VIRGO experiment 
[2]). These experiments should detect the gravitational waves generated by inspiralling 
compact binary systems at cosmological distances [3]-[7]. These are systems of two 
compact objects (neutron stars or black holes) spiralling very rapidly around each 
other in their late stage of evolution immediately preceding the final coalescence. 
The dynamics of such systems is entirely ruled by the radiation reaction forces due 



1 



to the emission of the gravitational radiation itself. Estimates of the number of 
inspir ailing/coalescing events each year have been found to be quite promising [8]- 
[11]. 

The LIGO and VIRGO observations of inspiralling compact binaries should 
provide precise measurements of the masses of the objects, possibly of their spins and 
possibly, in the case of neutron stars, of their radii [12]-[19]. The absolute luminosity 
distance of the binary will be measured independently of any assumption concerning 
the nature (and masses) of the objects [20]. But if the intrinsic masses of the objects 
are known, the cosmological red-shift of the host galaxy where the event took place 
can be measured. There is hope to deduce in this way a measurement of the Hubble 
parameter from an assumption concerning the statistical distribution of masses of 
neutron stars [21], [22]. The puzzle of the origin of gamma-ray bursts could be solved 
by comparison of the times of arrival of the gravitational waves and of gamma ray 
bursts [10], [23], [24]. Furthermore new limits on the validity of alternative theories of 
gravity, notably scalar-tensor theories [25], and new tests of general relativity in the 
strong field regime [26] should be possible by monitoring very precisely the inspiralling 
signal. 

To prepare for the analysis of such signals in the LIGO and VIRGO detectors 
one needs to compute the gravitational radiation field generated by a system of two 
point-masses moving on a circular orbit (the relevant case because the orbit will have 
been circularized by radiation reaction forces). Since inspiralling compact binaries 
are very relativistic, this problem is highly nontrivial and represents a challenge to 
relativity theorists [5], [7], as its resolution involves carrying out the calculation to a 
very high order in terms of a post- Newtonian expansion (see Will [27] for a review). 
The problem can be decomposed into two different problems, which can be referred to 
respectively as the "wave generation problem" and the "radiation reaction problem" 
(see Damour [28] for a discussion). 

The wave generation problem deals with the computation of the gravitational 
waveforms generated by the binary (at the leading order in 1/r, where r is the distance 
of the binary) when the orbital phase and frequency of the binary take some given 
values (j) and uj. This problem involves computing the amplitude of each harmonic 
of the wave corresponding to frequencies which are multiples of the orbital frequency, 
with the predominant harmonic being at twice the orbital frequency. 
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The radiation reaction problem consists of determining the evolution of the orbital 
phase (p(t) itself as a function of time, from which one deduces the orbital frequency 
u(t) = d(j)(t)/dt. The actual time variation of <f>(t) is nonlinear because the orbit 
evolves under the effects of gravitational radiation reaction forces. In principle it 
should be determined from the knowledge of the radiation reaction forces acting locally 
on the orbit. However these forces are at present not known with sufficient accuracy 
(only the relative first post- Newtonian corrections are known [29], [30]), so in practice 
the phase evolution is determined by equating the high-order post-Newtonian energy 
flux in the waves or energy loss (averaged over one orbit) and the decrease of the 
correspondingly accurate binding energy of the binary. 

Estimates of the precision needed in the resolution of these two problems can be 
inferred from black-hole perturbation techniques in the special case where the mass 
of one object is very small as compared with the other mass [31]-[37]. The required 
precision is reached when the systematic errors due to the neglect of some higher-order 
approximation become less than the statistical errors due to noise in the detector. It 
turns out that the post-Newtonian corrections in the time evolution of the phase 
(radiation reaction problem) will be measurable in advanced detectors, probably up 
to three post-Newtonian (3PN) orders beyond the (Newtonian) quadrupole radiation 
[33] , [37] . This corresponds to relativistic corrections in the energy loss as high as order 
~ (v/co) 6 where v is the orbital velocity and where we denote for later convenience 
the speed of light by cq. The possibility of measuring such high-order corrections can 
be understood crudely from the fact that, in order not to suffer a too severe reduction 
in signal-to-noise ratio, one will have to monitor the phase evolution with an accuracy 
of one tenth of a cycle over the tens of thousands of cycles during the entire passage 
through the frequency bandwidth of the detector. 

It has been argued [15]-[17] that most of the accessible information allowing 
accurate measurements of the binary's intrinsic parameters (such as the two masses) is 
in fact contained in the phase, because of the accumulation of cycles, and that rather 
less accurate information is available in the wave amplitude itself. For instance, the 
relative precision in the determination of the distance r of the source, which affects 
the wave amplitude, will be less than for the masses, which strongly affect the phase 
evolution. Consequently, the determination of r does not necessitate as high a post- 
Newtonian precision as for the masses, and to a good approximation it should be 
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adequate for that purpose to include in the gravitational waveform (as determined 
by solving the generation problem) post-Newtonian corrections up to an order less 
than 3PN order. To determine just what order is required in the waveform would 
necessitate a measurement-accuracy analysis similar to the ones performed in Refs. 
[17]-[19] but not restrained to a model waveform taking only the higher-order phase 
evolution into account (namely the so-called restricted post- Newtonian waveform). 
To our knowledge no complete investigation has been performed yet, however we 
believe that such a study would be desirable. For instance, it is conceivable that 
conclusions reached regarding the needed accuracy in the determination of the phase 
in the reaction problem are modified when the full amplitude evolution of the wave is 
taken into account. 

Both the gravitational waveform and associated outgoing flux or energy loss have 
been obtained recently to the 2PN approximation corresponding to the relative post- 
Newtonian order ~ (v/cq) 4 - Two computations have been performed, one by Blanchet, 
Damour and Iyer [38] using a mixed post-Newtonian and post-Minkowskian formalism 
[39], [40], and one by Will and Wiseman [41] using a purely post-Newtonian formalism 
due to Epstein and Wagoner [42] and Thorne [43]. The result of these computations 
for the energy loss has been summarized in [44]. More recently, the precision on the 
energy loss has been extended to include the next 2.5PN approximation [45]. 

Even though the 2PN or even 2.5PN precision in the resolution of the reaction 
problem appears to be still insufficient to make full use of the phase data [33], [37], 
it is plausible, as previously argued [15]-[17], that the 2PN waveform amplitude is 
already close to what will be needed by the LIGO and VIRGO detectors. However 
this 2PN waveform is displayed in Refs. [38], [41] in a format which is not ready 
for use in the future analysis of the outputs of LIGO and VIRGO. The aim of this 
paper is to present the complete, ready-to-use expression of the waveform, including 
the amplitude of all the wave harmonics present, to 2PN order, for non-spinning 
bodies. More precisely, we compute the two independent polarization states of the 
gravitational wave (customarily referred to as the "plus" and "cross" polarizations) 
which define the theoretical "templates" to be cross correlated with the raw output 
of the detectors. The templates should be used both "on line" when searching for the 
signal and later when the signal is subject to a very precise data analysis involving 
a more accurate determination of the parameters, and possibly the measurement of 
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other parameters. 

In Section 2 we present our main result, namely the two plus and cross polarization 
waveforms to 2PN order as functions of the phase and frequency. For completeness 
and discussion thereafter we recall also the results obtained in [38], [41] concerning 
the time evolution to 2PN order of the phase and frequency. In Section 3 we discuss 
various issues associated with the actual use of these waveforms. 

2. The gravitational waveforms to second-post-Newtonian order 

The gravitational wave as it propagates through the detector in the source's wave 
zone is entirely described by the transverse and traceless (TT) asymptotic waveform 
hJj T = (gij — $ij) TT where denotes the spatial covariant metric in a coordinate 
system adapted to the wave zone, and dij is the spatial part of the Minkowski 

metric (signature h ++)• The two polarization states h + and h x are defined 

by h + = \{piPj — qiqj)hj^ and h x = ^(piqj + Pjq^hJ? where p and q denote 
two polarization vectors forming, along with the direction n from the source to the 
detector, an orthonormal right-handed triad. The detector is directly sensitive to that 
linear combination of polarization waveforms h + and h x which is given by 

h(t) = F+h+(t) + F x h x (t) , (1) 

where F + and F x are the so-called beam-pattern functions of the detector depending 
on two angles giving the direction — n of the source as seen from the detector and 
a polarization angle specifying the orientation of the vectors p and q around that 
direction. The expressions of F + and F x in terms of these angles are given explicitly 
in the case of laser-interferometer detectors by e.g. Eq. (104) of Thorne [3]. 

The two polarizations h+ and h x in the case of a binary made of non-spinning 
point- masses moving on a (quasi-) circular orbit are obtained by a straightforward 
computation starting from the end results of Refs. [38], [41]. We choose the 
polarization vectors p and q to lie along the major and minor axis (respectively) 
of the projection onto the plane of the sky of the circular orbit, with p oriented 
toward the "ascending node" , the point at which body 1 crosses the plane of the sky 
moving toward the observer. The TT waveform hj? is split as in Eq. (4.1) of Ref. [38] 
into the sum of an "instantaneous" contribution (/i^ T )i ns t and of a "tail" contribution 
(/i^ T )taii (we comment below on the presence of wave tails in the signal and their 
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possible detection). The instantaneous contribution is given e.g. by Eqs. (4.5) of 
Ref. [38] from which we compute first (/i + ) inst and (/i x )inst- Then we add to these 
the tail parts (/i+, x )taii which are already obtained in Eqs. (4.9)-(4.10) of Ref. [38]. 
Identical formulae result from calculating h+ and h x from the circular-orbit limit of 
Eqs. (6.10) and (6.11) of Ref. [41]. Some transformations are necessary in order to 
express the result in the best form for applications in the LIGO and VIRGO detectors. 
The final result reads 

= ^ (^) 2/3 K> + *^<>? + *u« + ^Hf!? + } , 

(2) 

where the brackets involve a post- Newtonian expansion whose various post- Newtonian 
terms are given for the plus polarization by 
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and for the cross polarization by 
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The notation is as follows. The post-Newtonian expansion in (2) is parametrized by 
x = (Gmcu/co) 2 / 3 where u> is the 2PN-accurate orbital frequency of the circular orbit 
(a; = 2n/P where P is the orbital period) and m = mi + m2 is the total mass of the 
binary. [The expansion (2) is valid up to the neglect of 2.5PN terms of order 0(x 5 / 2 ).] 
In addition to m we denote dm = m± — m2 and rj = m\milm 2 . The vector n along 
the line of sight from the binary to the detector defines the inclination angle % with 
respect to the normal to the orbital plane. The normal is chosen to be right-handed 
with respect to the sense of motion so that we have < i < tt. The notations c and s 
are shorthand for the cosine and sine of the inclination angle : c = cosi and s = sinz. 



Finally the basic phase variable ip entering (3)- (4) is defined by 



if) = tj>- 



2Gmu 
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where <fi is the actual orbital phase of the binary, namely the angle oriented in the 
sense of motion between the ascending node and the direction of body 1(0 = mod 2n 
when the two bodies lie along p, with body 1 at the ascending node). The logarithmic 
term in the definition of ip involves a constant frequency scale uo which can be 
chosen arbitrarily (see below). This logarithmic phase modulation was determined 
in Refs. [46], [47] and originates physically from the propagation of tails in the wave 
zone. Using ift instead of the actual phase <p simplifies somewhat the expression of the 
waveform, since it permits collecting in a single block all the logarithmic terms (to 2PN 
order) . The variable ij) is also very convenient in black-hole perturbation theory, where 
it permits the resolution of the Teukolsky equation governing the outgoing radiation 
up to a very high order in the post- Newtonian expansion [35] , [36] . 

The expressions (2)-(5) solve the generation problem for inspiralling compact 
binaries to 2PN order. Up to now, both the plus and cross polarization waveforms 
were known for arbitrary masses to 1.5PN order [47], [48], and in the test mass limit 
7] -> to 1.5PN order [32] and 4PN order [36]. Through 1.5PN order equations (3) 
and (4) agree with Wiseman [49] (see also [27]) when those formulae are reduced 
to the circular-orbit limit. Equations (3) and (4) also largely agree with similar 
formulae given by Poisson [32] and Wiseman [47], however those formulae contained 
some typographical errors (see [48]). In any case, the final results given in Eqs. (3) 
and (4) supersede the previous results. We have checked that our expressions (3) 
and (4) in the test mass limit r\ — > are in perfect agreement with the truncation 
to 2PN order of the results given in Appendix B of Tagoshi and Sasaki [36]. (The 
comparison shows that the phase variable used by Tagoshi and Sasaki is related to ours 
by V'ts = ip + (2Gmo;/co)[ln2 — i|] which clearly corresponds simply to a rescaling of 
the frequency u>o.) 

In order to obtain the expressions of the waveforms h + (t) and h x (t) as functions 
of time one needs to replace ip and u appearing in equations (2)-(5) by their explicit 
time evolutions ip(t) and u(t) obtained from the resolution of the radiation reaction 
problem. This problem has been solved in Refs. [38], [41] to 2PN order, and we 
quote here these results for completeness and later discussion. It is convenient to 
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introduce instead of the local time t flowing in the experimenters' reference frame the 
dimensionless time variable 

e = ^fc-«), (6) 

where t c is a constant which represents the instant of coalescence of the two point- 
masses (at which the frequency goes formally to infinity). Then the instantaneous 
orbital phase <f>(t) is given in terms of the time variable (6) by Eq. (4.29) in Ref. [38] 
which reads 



, (() = , c -I{e^ + (^ + g,)e3/s-^ 



/ 9275495 284875 1855 2 \ 1/8 1 
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where <p c is another constant representing the value of the phase at t c . [Note that when 
taking into account higher order post-Newtonian approximations (starting at 2.5PN) 
the phase 4>(i) no longer tends to a constant when t — > t c but instead becomes infinite 
[45]. In this case the constant (fi c is simply determined by initial conditions when 
the frequency of the wave enters the detector's bandwidth.] The orbital frequency 
is obtained simply by differentiating equation (7) with respect to time (u> = d(f)/dt), 
hence 

( 1855099 56975 371 9 \ n 7/R ) 

+ H ??H V © 7 > ■ (8) 

V 14450688 258048' 2048 ' J J w 

From (7) and (8) one deduces the phase variable ip using equation (5). Both the 
expressions (7) and (8) are valid up to the 2PN order which corresponds formally 
to the same relative precision as for the waveforms (3) and (4). However it is not a 
priori required for consistency to have the same post- Newtonian precision in both the 
waveforms and phase. On the contrary, one should use in the waveforms (3) and (4) 
the best available expression for <p(t) which will be hopefully determined in the future 
to a much higher order than 2PN (see Ref. [45] for the expression of <p(t) to 2.5PN 
order). Related to this, note that the logarithmic term in the phase variable t(j(t) 
although of formal order 1.5PN, is actually of order 4PN relatively to the dominant 
term in (f)(t) given by (7) (indeed 4>{t) — 4> c is of order ~ Cq which is the inverse of the 
order of the radiation reaction effects) . Thus the logarithmic term is in fact currently 
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negligible but must be included when the precision on <f>(t) improves in the future to 
reach 4PN (even without knowing <f>(t) to the 4PN order it may be a good idea to 
include this term in the templates). 

It can be readily checked from the expressions (5)-(8) that any change in the 
phase variable ip corresponding to a rescaling of the frequency c^o is equivalent, to 
the considered order, to a shift in the value of the instant of coalescence t c , namely, 
the rescaling uj g — > Acu is equivalent to the shift t c — > t c — (2Gm/co)lnA. Thus 
a particular choice of the frequency scale ujo is physically irrelevant since a different 
choice can always be absorbed into a redefinition of the origin of time in the wave zone. 
(We have u = (l/46)exp(^| — C) where C = 0.577... is Euler's constant and b is a 
freely specifiable parameter entering the relation between the wave-zone coordinate 
time t and the harmonic coordinates tn,rn '■ t = tjj — th/cq — (2Gm/c^) ln(r#/co&) 
where th is the distance of the source in harmonic coordinates, see [46].) It could be 
possible to relate ujq to the source characteristics, choosing for instance ujq = Cq/Gui 
where m is the total mass of the binary. However this mass will not be known in 
practice but will be used as a parameter in the templates to be varied during the data 
analysis process, so this choice seems to be somewhat awkward. For practical purposes 
it is probably more convenient to choose too in such a way that it is uniform over all 
templates. For instance one can relate u to the detector characteristics by choosing 
ujo/tv = 10 Hz where 10 Hz (say) is the seismic cutoff frequency of the detector [26]. 
Here we adopt this choice. 

3. Discussion 

To construct adequate filters for the analysis of inspiralling binary signals, one 
should proceed as follows. The theoretical waveform h(t) given as a function of 
time by the above formulas (l)-(8) is discretized and its Fourier transform h(Q) 
is computed numerically and stored. Then the ratio q(ft) = h(fl)/S n (ft), where 
S n (Cl) is the measured power spectral density of the noise in the detector, defines the 
Fourier transform of the Wiener filter q(t), which is to be cross- correlated (either in 
real time or during the more precise analysis later) with the raw output o(t) of the 
detector composed of the superposition of the actual signal and of the noise (which is 
supposed here to be Gaussian). Because of the availability of fast Fourier transforms 
the correlation is computed in the Fourier domain using the (discrete) correlation 

10 



theorem for the Fourier coefficients q(fl) and o(Q). This is repeated for each set of 
parameters in the filter until maximization of the signal-to-noise ratio is obtained 
yielding the determination of the parameters of the binary (if a signal was really 
present at this instant). 

For two nonrotating test-masses there are 4 parameters entering the phase of the 
signal (this is true up to any post-Newtonian order). These can be chosen to be the 
two mass parameters m and rj, and the constants t c and (j) c . Alternatively one can use 
the chirp mass = rf^m and the reduced mass \i = rim, and/or the arrival time to 
and phase 4>o at the seismic frequency u$. The mass parameters may also be replaced 
by two distinct post-Newtonian pieces of the (chirp) time left till coalescence starting 
from ujq [50]. In addition, the amplitude of each harmonic of the signal depends on 
the distance r, on the inclination angle i, and on the direction of the binary and the 
polarization angle through the beam-pattern functions F + and F x . Note that r is the 
cosmological luminosity distance, and that the masses are the red-shifted masses which 
are equal to the true masses multiplied by 1 + z where z is the binary's cosmological 
redshift. This permits the cosmological measurements proposed in Refs. [20]-[22]. 

In the case of rotating bodies, there are additional parameters in the signal. Spin- 
orbit and spin-spin contributions to the waveform and the phase evolution have been 
obtained by Kidder et al [51], [52]. The over-all effect of spins on the accumulated 
phase during an inspiral was analyzed by Blanchet et al [44]. This calculation was 
done with the assumption that the spins of the objects were aligned and perpendicular 
to the orbital plane. However, including non-aligned spins makes the waveform and 
accumulated-phase calculations considerably more complicated. In this general case 
the orbital plane can wobble - thus the inclination angle i in Eqs. (3) and (4) changes 
with time - giving rise to an amplitude modulation and a frequency modulation [53]. 
However note that for realistic inspiralling neutron star binaries (such as the binary 
pulsar PSR 1913+16 at the epoch when it finally coalesces) the spin-orbit and spin- 
spin parameters are expected to be very small, so that the dynamics of these systems 
will be dominated by the purely gravitational effects investigated here (see discussion 
in [44]). 

In addition to involving rotating bodies, the binary could move on a slightly 
eccentric orbit if, being formed late, it reaches the final inspiral stage before the 
circularization of the orbit by radiation reaction has fully taken place. This would 
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introduce an additional term in the phase (7) involving a new parameter which can 
conveniently [19] be taken to be Cqo;q 9//9 where eo is the initial eccentricity at uo (recall 
that e 2 evolves in time proportionaly to cu _19//9 in the quadrupole approximation). 
The general formulae of Refs. [40], [41] and [54] can be applied to such non-circular 
situations. 

Another possibility is that general relativity is not the correct theory of gravity 
because, for instance, of the existence of a scalar spin-0 field besides the spin-2 metric 
field. In this case there would be still another parameter in the signal describing the 
coupling of the scalar field with the matter fields. Preliminary analyses of the bounds 
that could be placed on the coupling parameter o;bd of the Brans-Dicke theory can 
be found in Refs. [25], [19]. 

More generally, it should be possible to check that the theoretical waveform (l)-(8) 
is exactly reproduced in the real signal without prejudice of which theory could be 
the true theory of gravity. A simple method to do this consists of parametrizing 
the templates by a redundant set of parameters (i.e. by more parameters than 
strictly necessary) and measuring them independently by optimal filtering. In this 
way one can test whether some particular terms involving some specific combinations 
of parameters are really present in the signal [26]. However since multiplying the 
number of independent parameters has the effect of diminishing the accuracy in their 
measurement, this method necessitates a high signal-to-noise ratio. 

An application of this method is the detection in the real signal of effects associated 
with the tails of waves. Physically the tails come from the backscattering of the linear 
waves off the static spacetime curvature generated by the total mass m of the binary 
itself. The tails are characterized by the fact that they are non "instantaneous", 
namely they depend on the whole integrated past of the source ; however for binary 
systems the actual dependence on the past is negligible to 2PN order. Nevertheless 
some important signatures of the tail effect remain in the signal, the most important 
of which is the term with coefficient 7r in the phase of the wave (7), which is directly 
due to the tails in the radiation reaction forces. This effect should be easily detectable 
in the real signal. Other tail contributions enter the binary's waveform and are due 
to the propagation of tails in the wave zone. An instance is the logarithmic phase 
modulation in (5) which can be referred to as a wave tail contribution to the phase 
(as opposed to the radiation reaction contribution). Even this very small tail-induced 
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phase modulation (of relative 4PN order) has been shown to be detectable with a 
high, but not exceedingly high, signal-to-noise ratio using the above method [26]. 
This shows how the observations of inspiralling compact binaries could permit new 
tests of general relativity in a regime of strong and rapidly varying gravitational fields. 
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